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By Ivan Nourdin and Giovanni Peccati 

Universite Paris VI and Universite Paris Ouest 

We show how to detect optimal Berry-Esseen bounds in the nor- 
mal approximation of functionals of Gaussian fields. Our techniques 
are based on a combination of Malliavin calculus, Stein's method and 
the method of moments and cumulants, and provide de facto local 
(one-term) Edgeworth expansions. The findings of the present paper 
represent a further refinement of the main results proven in Nour- 
din and Peccati [Probab. Theory Related Fields 145 (2009) 75-118]. 
Among several examples, we discuss three crucial applications: (i) 
to Toeplitz quadratic functionals of continuous-time stationary pro- 
cesses (extending results by Ginovyan [Probab. Theory Related Fields 
100 (1994) 395-406] and Ginovyan and Sahakyan [Probab. Theory Re- 
lated Fields 138 (2007) 551-579]); (ii) to "exploding" quadratic func- 
tionals of a Brownian sheet; and (iii) to a continuous-time version 
of the Breuer-Major CLT for functionals of a fractional Brownian 
motion. 

1. Introduction. Let {F n : n > 1} be a sequence of zero-mean real-valued 
random variables and consider a standard Gaussian variable N ~ ^V(0, 1). 
Assume that each F n is a functional of an infinite-dimensional Gaussian field 
and suppose that, as n — > oo, 

(1.1) F n ^N. 

In the paper [20], the present authors demonstrated that one can naturally 
combine Malliavin calculus (see, e.g., [13, 21]) with Stein's method (see, e.g., 
[4, 29, 33, 34]) in order to obtain explicit bounds of the type 

(1.2) d(F n ,N)<ip(n), n>l, 
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where d(F n ,N) stands for some appropriate distance (e.g., the Kolmogorov 
distance or the total variation distance) between the laws of F n and N, and 
<p(n) is some positive sequence converging to zero. The aim of the present 
work is to develop several techniques, allowing us to assess the optimality of 
the bound </?( n ) appearing in (1.2) for a given sequence {F n }. Formally, one 
says that the bound </?(n) is optimal for the sequence {F n } and the distance 
d whenever there exists a constant c G (0, 1) (independent of n) such that, 
for n sufficiently large, 

(1.3) c<d(F n ,N)/<p(n)<l. 

We shall establish relations such as (1.3) by pushing the Malliavin-type ap- 
proach to Stein's method (initiated in [20]) one step further. In particular, 
the findings of this paper represent a new and substantial refinement of the 
central limit theorems (CLTs) for functionals of Gaussian fields which were 
proven in [22, 23, 25, 26]. Once again, our techniques do not require that the 
random variables {F n } have the specific form of partial sums. Indeed, we 
will see, in Sections 4-6 below, that our results yield optimal Berry— Esseen- 
type bounds for CLTs involving objects as diverse as: (i) Toeplitz quadratic 
functionals of continuous-time stationary processes; (ii) quadratic function- 
als of a Brownian motion or of a Brownian sheet indexed by a compact set 
of M. d (d>2); and (hi) polynomial functionals constructed from a fractional 
Brownian motion. 

Note that, in the subsequent sections, we shall focus uniquely on the 
normal approximation of random variables with respect to the Kolmogorov 
distance. This distance is defined as 

(1.4) d Kol (X, Y) = sup\P(X <z)- P(Y < z)\ 

for any pair of random variables X and Y . It will later become clear that 
many results of the present paper extend almost verbatim to alternate dis- 
tances, such as the Wasserstein and the total variation distance, between 
laws of real- valued random variables. 

Our basic approach can be described as follows. Fix z € M. and consider 
the Stein equation 

(1.5) l { _ 00jZ] (x)-<f>(z)=f(x)-xf(x), x£R, 

where, here, and for the rest of the paper, we use the standard notation 
$(z) = P(N < z) [N ~ ^Y(0, 1)] and where 1a stands for the indicator of a 
set A. It is well known that, for every fixed z, (1.5) admits a solution f z such 
that ||/z||oo <■ \/27r/4 and H/^Hoo < 1 (see, e.g., [4], Lemma 2.2, or formulae 
(2.10)-(2.11) below). Now, suppose that the elements of the sequence {F n } 
appearing in (1.1) are functionals of some Gaussian field, say X, and assume 
that each F n is differentiable in the sense of Malliavin calculus (see Section 
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2.1 for details). Denote by DF n the Malliavin derivative of F n and by L^ 1 the 
pseudo-inverse of the Ornstein-Uhlenbeck generator (again, see Section 2.1). 
Recall that DF n is a random element with values in an appropriate Hilbert 
space fj. In [20], Section 3, we proved and applied the following relations, 
direct consequences of the fact that f z solves (1.5) on the one hand, and 
of the celebrated integration by parts formula of Malliavin calculus on the 
other hand: for every zGl, 

P(F n < z) - <f>(z) = E[f' z (F n ) - F n f z (F n )] 

(1.6) 

= E[f' z (F n )(l - (DF n , -DL^FM]. 

By using (1.4), applying the Cauchy-Schwarz inequality to the right-hand 
side of (1.6) and using the fact that f' z is bounded by 1, one immediately 
obtains that 



(1.7) d Kol {F n ,N)<^E[{l-{DF n ,-DL-^F n )^f]. 

The starting point of [20] was that, in several crucial cases (e.g., when each 
F n is a multiple Wiener-Ito integral of a fixed order), the upper bound 



(1.8) p(n):=^E[(l-(DF n ,-DL-iF n )z)2], n>l, 

is such that: (i) the quantity ip(n) can be explicitly computed (e.g., in terms 
of contraction operators); (ii) (f(n) — ► as n — ► oo; and (iii) <p(n) is directly 
related to quantities playing a fundamental role in the CLTs for functionals 
of Gaussian fields proven in [22, 23, 25, 26]. The aim of the present paper is 
to establish conditions on the sequence {F n } ensuring that the ratios 

(1.9) WUF.M-VF-.-DL-iFM] 

(p[n) 

involving (1.8) and the right-hand side of (1.6), converge to a nonzero limit 
for all z outside some finite set. Such a result immediately yields the exis- 
tence of a constant c, verifying (1.3) for d = d^ \. We will show that a very 
effective way to prove the convergence of the quantities appearing in (1.9) is 
to characterize the joint convergence in distribution of the random vectors 

(1.10) [F n , j, n>l, 

toward a two-dimensional Gaussian vector with nonzero covariance. The ap- 
plications presented in Sections 4-6 will show that this specific convergence 
takes place in several crucial situations, involving, for instance, quadratic or 
polynomial functionals of stationary Gaussian processes. We will see that, in 
order to prove a CLT for the vector appearing in (1.10), a useful tool is the 
multidimensional version of the CLT for multiple stochastic integrals which 
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was proven in [26]. Also, it is interesting to note that if each F n in (1.1) is a 
double stochastic integral, then our conditions can be expressed exclusively 
in terms of the second, third, fourth and eighth cumulants associated with 
the sequence {F n }; see Section 3.3 below. 

The rest of the paper is organized as follows. Section 2 deals with pre- 
liminaries concerning Malliavin calculus, Stein's method and related topics. 
Section 3 contains our main results, with special attention devoted to ran- 
dom variables belonging to the second Wiener chaos of a Gaussian field. 
In Section 4, we develop an application to Toeplitz quadratic functional of 
stationary continuous-time Gaussian processes, thus extending and refining 
some results by Ginovyan [7] and Ginovyan and Sahakyan [8]. Section 5 is 
devoted to quadratic functionals of Brownian motion and of the Brownian 
sheet, whereas Section 6 focuses on a continuous-time version of the Breuer- 
Major CLT for processes subordinated to a fractional Brownian motion. 

2. Preliminaries. 

2.1. Gaussian fields and Malliavin calculus. We shall now provide a 
short description of the tools of Malliavin calculus that will be needed in 
the forthcoming sections. The reader is referred to the monographs [13] and 
[21] for any unexplained concepts or results. 

Let 9) be a real separable Hilbert space. We denote by X = {X(h) : h £ 
fj} an isonormal Gaussian process over Sj. By definition, X is a centered 
Gaussian family indexed by the elements of f) and such that, for every 
h,gefi, 

(2.1) E[X(h)X(g)] = (h,g) Si . 

In what follows, we shall use the notation L 2 (X) = L 2 (f2, a(X),P). For every 
q > 1, we write 9)® q to indicate the qth tensor power of Sj; the symbol Sj <3q 
stands for the qth symmetric tensor power of S), equipped with the norm 
y/ql\\ • We denote by I q the isometry between 9j Gq and the gth Wiener 
chaos of X. It is well known (again, see [21], Chapter 1, or [13]) that any 
random variable F belonging to L 2 (X) admits the chaotic expansion 

oo 

(2-2) F = J2W q ), 

g=0 

where Io(/o) := -^[-^li the series converges in L 2 and the kernels f q E S) Qq , 
q > 1, are uniquely determined by F. In the particular case where 9) = 
L 2 (j4, fj,), where (A,g/) is a measurable space and fi is a <7-finite and 
nonatomic measure, one has that Sj Qq = L 2 (^4 9 , £/® q , n® q ) is the space of 
symmetric and square-integrable functions on A q . Moreover, for every / £ 
Sj &q , I q (f) coincides with the multiple Wiener-Ito integral (of order q) of 
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/ with respect to X (see [21], Chapter 1). It is well known that a random 
variable of the type I q (f), f G $) &q , has finite moments of all orders (see, e.g., 
[13], Chapter VI). Moreover, any nonzero finite sum of multiple stochastic 
integrals has a law which is absolutely continuous with respect to Lebesgue 
measure (see, e.g., Shigekawa [32] for a proof of this fact; see [21], Chapter 1, 
or [30] for a connection between multiple Wiener-Ito integrals and Hermite 
polynomials on the real line). For every q > 0, we denote by J q the orthogonal 
projection operator on the qth. Wiener chaos associated with X so that, if 
F G L 2 (<rP0) is as in (2.2), then J q F = I q (f q ) for every q > 0. 

Let {efc,fc > 1} be a complete orthonormal system in Sj. Given / G fj 0p 
and g G S) &q , for every r = 0, . . . ,p A q, the rth contraction of / and g is the 
element of f)®(p+i- 2r ) defined as 

oo 

(2.3) f® r 9= ^' e *i ® "' ® Oft®- ® {9,e h ® ••• O^r. 

il,...,i r =l 

In the particular case where $j = L 2 (A,£/,/j,) (with \i nonatomic), one has 
that 

f® r g= / f(ti,...,tp- r ,Si,...,S r ) 

x g(tp-r+l, tp+q-2r, Si,..., S r ) dfl(si) ■ ■ ■ d/J,(s r ). 

Moreover, / ®o 9 = / ® 9 equals the tensor product of / and g while, for 
p = q, f ® p g = (f,g)fi®p. Note that, in general (and except for trivial cases), 
the contraction f® r g is not a symmetric element oi 9)® <yP+q ~ 2r \ The canon- 
ical symmetrization of f ® r 9 is written f<S> r g. We also have the following 
multiplication formula: if / G f)® p and g G S)® q , then 

pAq f \ / \ 

(2.4) Ip{f)I q (9) =5> ! W2r(/M- 

Let ^ be the set of all smooth cylindrical random variables of the form 

F = g(X{ ( f> l ),...,X(<p n )), 

where n > 1 , g : M n — ► R is a smooth function with compact support and 
<j>% G Sj. The Malliavin derivative of F with respect to X is the element of 
L 2 (0,fj) defined as 

i=l * 

Also, DX(h) = h for every ft G f). By iteration, one can define the mth 
derivative D m F [which is an element of L 2 (Q,$)® m )] for every m>2. As 
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usual, for m > 1, B m ' 2 denotes the closure of 5? with respect to the norm 
II • \\m,2, defined by the relation 

m 

1=1 

Note that if F is equal to a finite sum of multiple Wiener-Ito integrals, then 
F G B" 1 ' 2 for every m > 1. The Malliavin derivative D verifies the following 
chain rule: if ip : W 2 — > K is in ^ (he., the collection of bounded continuously 
differentiable functions with a bounded derivative) and if {i*i}j=i n is a 
vector of elements of B 1,2 , then ip(Fi, . . . , .F n ) G B 1,2 and 

D<p(F 1 ,...,F n )=J2-^ : (F 1 ,...,F n )DF i . 

i=l 1 

Observe that the previous formula still holds when ip is a Lipschitz func- 
tion and the law of (F\, . . . , F n ) has a density with respect to the Lebesgue 
measure on W 1 (see, e.g., Proposition 1.2.3 in [21]). We denote by 5 the 
adjoint of the operator D, also called the divergence operator. A random 
element u G L 2 (0,ij) belongs to the domain of 5, noted Dom.5, if and only 
if it verifies 

\E{DF,u)^\<c u \\F\\v for any Fe/, 

where c u is a constant depending uniquely on u. If u G Dom<5, then the ran- 
dom variable <5(u) is defined by the duality relationship (i.e., the "integration 
by parts formula") 

(2.5) E(F8{u))=E(DF,u)%, 

which holds for every FgD 1 ' 2 . 

The operator L, acting on square-integrable random variables of the type 
(2.2), is defined through the projection operators {J q } q >o as L = J2q^=o ~Q.Jq 
and is called the infinitesimal generator of the Ornstein-Uhlenbeck semi- 
group. It verifies the following crucial property: a random variable F is an 
element of DomL (=B 2,2 ) if and only if F G Bom5D (i.e., F G B 1,2 and 
DF G Dom<5) and, in this case, 5DF = —LF. Note that a random variable 
F as in (2.2) is in B 1 ' 2 (resp., B 2 ' 2 ) if and only if 

oo / oo \ 

^<l\\fq\\%Q q < oo resp., J2q 2 \\fq\\l@ q < oo , 

q=l \ q=l I 

and also E^H-D-FH^] = J2q>i l\\fq\\%®q- If # = L 2 (A, £/ , fi) (with \i nonatomic), 
then the derivative of a random variable F as in (2.2) can be identified with 
the element of L 2 (A x Q) given by 

oo 

(2.6) D a F = Y^qI q -i{f q (;a)), at A. 

q=l 
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We also define the operator L" 1 , which is the pseudo-inverse of L, as follows: 
for every F G L 2 (X), we set L~ X F = Y. q >i\ J q( F )- Note that L_1 is an 
operator with values in D 2 ' 2 and that LL~ l F = F — E(F) for all F G L 2 (X). 
The following lemma generalizes Lemma 2.1 in [19]. 

Lemma 2.1. Xei F G B 1 ' 2 be such that E(F) = 0. Suppose that, for some 
integer s>0, E\F\ s+2 < oo. Then, 

(2.7) E(F S (DF, -DL- 1 ^) = j^-E(F s+2 ). 



PROOF. Since L 1 F G P 2 ' 2 , we can write 
-DL^F)^ 

~E({D(F s+1 ),D(-L~ 1 F))f ) ) 



E(F S (DF, —DL^F)^) 
1 



□ 



a + 1 

= l —E{5DL~ 1 F x F s+1 ) [by integration by parts (2.5)] 

s ~\- \ 

= -^—E(F S+2 ) (by the relation —8DL~ X F = F). 
s + 1 

Remark 2.2. If F = I q (f) for some q>2 and / G £> 09 , then 

(2.8) {DF-DL- l F)„ = (DI q (f), -DL~ l I q {f)) ^ = -\\DI q (f)\\ 2 

so that (2.7) yields, for every integer s > 1, that 

(2.9) S(J g (/)'||D/,(/)||J) = ^- l E(I g (fy+ 2 ). 

2.2. Stein's method and normal approximation on a Gaussian space. We 
start by recalling that, for every fixed z G K, the function 

(2.10) h{x) = e* 2 / 2 f [l^^a)-*^^^ 

J — oo 

/ 211 n = f^/2¥e a;2 / 2 $(x)(l-$( Z )), ifa<z, 

1 ' ' l^ 2 / 2 $(z)(l-$(i)), ifx>z, 



is a solution to the Stein equation (1.5), also verifying ||/ 2 ||oo < v27r/4 and 

H/illoo < 1- 

The following lemma will play a crucial role in the sequel; see also (1.6). 
Its content is the starting point of [20]. 
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Lemma 2.3. Let F E B 1,2 have zero mean. Assume, moreover, that F 
has an absolutely continuous law with respect to the Lebesgue measure. Then, 
for every z E R, 

P(F <z)- = E[f' z (F)(l - (DF, -DL^F)*)}. 

Proof. Fix zeK. Since f z solves the Stein equation (1.5), we have 
P(F < z) - $0) = E[f' z (F) - Ff z (F)]. Now, observe that one can write 
F = LL~ 1 F = —5DL~~ 1 F. By using the integration by parts formula (2.5) 
and the fact that Df z (F) = f' z {F)DF (note that, for this formula to hold 
with f z only Lipschitz, one needs F to have an absolutely continuous law — 
see Section 2.1), we deduce that 

E[Ff z (F)}=E[-5DL- l Ff z (F)} 

= E[{Df z (F),-DL~ 1 F)d 

= E[f' z (F)(DF,-DL- 1 F)^}. 

It follows that E[f' z (F) - Ff z (F)] = E[f' z (F)(l - {DF, -DL^F)^)} and the 
proof of the lemma is complete. □ 

As an application, we deduce the following result, first proven in [20] (the 
proof is reproduced here for the sake of completeness). 

Theorem 2.4. Let F e B 1,2 have zero mean and N ~ ^(0, 1). Then, 

(2.12) d Kol (F, N) < ^E[{l-{DF,-DL-^F)„Y]. 

IfF = L q (f) for some q>2 and f G Sj Q i , then (DF, -DL~ X F)^ = q-^DF^ 
and therefore 

(2.13) d Kol (F, N) < ^/ E[(l - q-i\\DF\\l) 2 }. 

Proof. If / is a bounded, continuously differentiable function such that 
ll/'lloo < 1) then, using the same arguments as in the proof of Lemma 2.3 
(here, since / belongs to observe that we do not need to assume that 
the law of F is absolutely continuous), we have 

\E[f'(F) - Ff(F)}\ = \E[f'(F)(l - (DF,-DL- l F)„)]\ 

< E\l - (DF,-DL^ l F)^\. 

In fact, the inequality \E[f'(F) - Ff(F)}\ < E\l - (DF.-DL^F)^ con- 
tinues to hold with f = f z (which is bounded and Lipschitz, with Lipschitz 
constant less than one) as is easily seen by convoluting f z by an approxi- 
mation of the identity. Hence, Lemma 2.3, combined with Cauchy-Schwarz 
inequality, implies the desired conclusion. □ 
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Remark 2.5. In general, the bound appearing on the right-hand side 
of (2.12) may be infinite. Indeed, the fact that F £ B 1,2 only implies that 
(DF,—DL~ 1 F)f ) € L 1 (J7). By using the Cauchy-Schwarz inequality twice, 
one sees that a sufficient condition, in order to have (DF, —DL~ 1 F)f ) £ 
L 2 (ft), is that \\DF\lf, and ||DL _1 F|| fi belong to L 4 (0). Also, note that 
if F is equal to a finite sum of multiple integrals (e.g., F is a polynomial 
functional of X), then the random variable {DF, —DL~ 1 F)f ) is also a finite 
sum of multiple integrals and therefore has finite moments of all orders. In 
particular, for F = Iq(f), the right-hand side of (2.13) is always finite. 

The bounds appearing in Theorem 2.4 should be compared with the forth- 
coming Theorem 2.6, dealing with CLTs on a single Wiener chaos (part A) 
and on a fixed sum of Wiener chaoses (part B). 

Theorem 2.6 (See [22, 23, 25, 26]). Fix q>2 and let the sequence 
F n = Iq{fn), n>l, where {f n } C Sj Gq , be such that E{F%] — ► 1 as n — ► oo. 

(A) The following four conditions are equivalent as n — > oo : 

(i) F n ^iV~^(0,l); 

(ii) E(F*)^3; 

(iii) \\f n <S)j f n \\fim<i-3) -> f° r ever V 3 = !> ■ • ■ ,1 ~ !/ 

(iv) l-q- l \\DF n \\l -+0 mL 2 . 

(B) Assume that any one of conditions (i)-(iv) of part A is satisfied. Let 
the sequence G n , n > 1, have the form 

M 

P =i 

for some M > 1 (independent of n) and some kernels g^t 1 £ Sj &p (p = 1, . . . , M, 
n > 1). Suppose that, as n — > oo, 

M 

£(G 2 J = ]>>!||^||^ P ^ C 2 >0 and \\g® ®j — 

P =i 

/or every p = 1, . . . , M and every j = 1, . . . ,p — 1. 7f i/te sequence of covari- 
ances E(F n G n ) converges to a finite limit, say i/ien (F n ,G n ) con- 

verges in distribution to a two-dimensional Gaussian vector (N\,N2) such 
that E{Nl) = I, E(N%) = c 2 and E{N 1 N 2 ) = p. 

The equivalence between points (i)-(iii) in part A of the previous state- 
ment was first proven in [23] by means of stochastic calculus techniques; the 
fact that condition (iv) is also necessary and sufficient was proven in [22]. 
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Part B (whose proof is straightforward and therefore omitted) is a conse- 
quence of the main results established in [25, 26]. Note that in part B of 

(p) 

the previous statement, we may allow some of the kernels g)i to be equal 
to zero. See [19] and [20], Section 3.3, for some extensions of Theorems 2.4 
and 2.6 to the framework of noncentral limit theorems. 

Remark on notation. In what follows, given two numerical sequences 
{a n } and {b n }, the symbol a n ~ b n means that lima n /6 n = 1, whereas a n x 
b n means that the ratio a n /b n converges to a nonzero finite limit. 

2.3. A useful computation. We shall denote by {H q : q > 0} the class of 
Hermite polynomials, defined as follows: Hq = 1 and, for q > 1, 

(2.14) H q (z) = {-iye z2 ' 2 ^ e~ z2 l\ z G R; 

for instance, H\(z) = z, H2(z) = z 2 — 1 and so on. Note that the definition 
of the class {H q } immediately implies the recurrence relation 

(2.15) A. Hq {z)e- z2 t 2 = -H q+1 (z)e~* 2 / 2 , 

yielding that the Hermite polynomials are related to the derivatives of 
= P{N <z) [JV~ ^(0, 1)], written $(«)(z) (q= 1,2,. . .), by the for- 
mula 



(2.16) $W(z) = (-iy- 1 H q _ 1 ( 



z)- 



-z 2 /2 



/2vr 

We also have, for any q > 1, 

d 

dz 



(2.17) —H q (z)=qH q _ 1 (z). 



Now, denote by f z the solution to the Stein equation (1.5) given in for- 
mulae (2.10)-(2.11). The following result, connecting f z with the Hermite 
polynomials and the derivatives of $, will be used in Section 3. 

PROPOSITION 2.7. For every q>l and every z £M., 

»+oo e ~ x2/2 1 e~ z2/2 



(2.18) 



L f '* ix)H « {x) ^ dX= —2 H « +l{z) -V2; 

= J_(_l)9+l$(9+2)( z ). 
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Proof. By integrating by parts and by exploiting relations (2.11) and 
(2.15), one obtains that 

r+oo e~~ x2 l 2 

/ f'z(x)H q (x) dx 

J — oo V Z7T 

(2.19) = / f z (x)H q+1 (x)^=^dx 

J—oo V Z7T 

/'-(-OO / px \ 

= —J^ H q+1 {x)^J^Jl { ^ z] {a)-^(z))e- a2 l 2 dayx. 

By integrating by parts, using H q+ \ = ^2^+2 [ see (2-17)] and in view of 
(2.15), one easily proves that 

f + °°H q+1 (x)( f (l(_ 00j , ] (a) -$>{z))e- a2 l 2 da\dx 

J—oo \J — 00 J 

1 r+oo „ 

= -— - / H q+2 {x){l { ^ z] (x)-${z))e-* > 2 dx 
Q ~r £ J—oo 



q + 
1 



J-( f if g+2 (x) e ^ 2 / 2 dx-$(z) [ +C ° H q+2 (x)e- x2 ' 2 dx) 

T £ \J —00 J — oo / 



-H q+ i{z)e 



z 2 /2 



q + 2 

By plugging this expression into (2.19), we immediately arrive at the desired 
conclusion. □ 



For instance, by specializing formula (2.18) to the case q = l, one obtains, 
for iV~^(0,l), 

(2.20) E[f>(N) x N] = \(z 2 - l)^f = i$(3) (z) . 

3. Main results. 



3.1. Two general statements. We start by studying the case of a general 
sequence of Malliavin derivable functionals. 

Theorem 3.1. Let F n , n > 1, be a sequence of centered and square- 
integrable functionals of some isonormal Gaussian process X = {X(h) :h G 
such that E(F 2 ) — ► 1 as n— >oo. Suppose that the following three con- 
ditions hold: 

(i) for every n, one has that F n £ D 1 ' 2 and F n has an absolutely contin- 
uous law (with respect to the Lebesgue measure); 
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(ii) the quantity ip{n) = \/E[(l - (DF n , - D L" 1 F n ) <^) 2 ] [as in (1.8)] is 
such that: (a) ip(n) is finite for every n; (b) as n — ► oo, ip{n) converges to 
zero; and (c) there exists m>l such that ip(n) > for n > m; 

(iii) as n — > oo, the two-dimensional vector (F n , ~~^~^7^ — ) / as 
in formula (1.10)] converges in distribution to a centered two-dimensional 
Gaussian vector (Ni,N 2 ) such that E(N?) = E(N%) = 1 and E(N 1 N 2 ) = p. 

Then, the upper bound d^ \(F n , N) < <p(n) holds. Moreover, for every zGK, 

(3.1) ^(n)- 1 ^ <z)- $(*)] _+J(^ - xfl^ = £**\z). 

As a consequence, if p^O, then there exists a constant c £ (0, 1), as well as 
an integer uq > 1, such that relation (1.3) holds for d = dxoi an d f or every 
n > no. 

Proof. Fix z € R. From assumption (i) and Lemma 2.3, recall that 

V {ny l [P{F n <z)- $>{z)]=E[f z (F n ) V {n)-\l - (DF n , -DL^F^)}. 

The facts that f' z is bounded by 1 on the one hand and that 9j(n) _1 (l — 
(DF n ,—DL~ 1 F n )f ) ) has variance 1 on the other hand imply that the se- 
quence 

fUFMn)- 1 ^ - (DF^-DL-'E^), n > 1, 

is uniformly integrable. Now, deduce from (2.10) that x — > f' z (x) is contin- 
uous at every x ^ z. This yields that, as n — > oo and due to assumption 
(iii), 

E\J' z {F n )ip(n)-\l - (DF n , -DL- l F n )^] 
— E(f' x (N\)N2) =px Etf'iN^N,). 

Relation (3.1) now follows from formula (2.20). If, in addition, p ^ 0, then 
one can obtain the lower bound (1.3) by using the elementary relation 

\P(F n <0)-H0)\<d Kol (F ni N). □ 

Remark 3.2. Clearly, the conclusion of Theorem 3.1 still holds when n 
is replaced by some continuous parameter. The same remark holds for the 
forthcoming results of this section. 

The next proposition connects our results with one-term Edgeworth ex- 
pansions. Note that, in the following statement, we assume that E{F n ) = 
and E(F%) = 1 so that the first term in the (formal) Edgeworth expansion of 
P{F n < z) — <3?(z) coincides with —(3\)~ 1 E(F^)^ 3 \z). For an introduction 
to Edgeworth expansions, the reader is referred, for example, to McCullagh 
[18], Chapter 3, or Hall [12], Chapter 2. See also Rotar [31] for another 
application of Stein's method to Edgeworth expansions. 
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Proposition 3.3 (One-term Edgeworth expansions). Let F n , n > 1, be 
a sequence of centered and square-integrable junctionals of the isonormal 
Gaussian process X = {X(h):h G $)} such that E{F%) = 1. Suppose that 
conditions (i)-(iii) of Theorem 3. 1 are satisfied and also that: 

(a) for every n, E\F n \ 3 < oo; 

(b) there exists e > such that sup n>1 E\F n \ 2+£ < oo. 
Then, as n — > oo, 

(3.2) _L_ £( ^>__ p 

and, for every z € K, one has the one-term local Edgeworth expansion 
(3.3) 



P(F n <z)- *(z) + ^E(F 3 )& 3 \z) = o s (<p(n)), 



where o z {ip{n)) indicates a numerical sequence (depending on z) such that 
V?(n) _1 o z (ip(n)) — > as n — > oo. 

Remark 3.4. Of course, relation (3.3) is interesting only when p ^ 0. 
Indeed, in this case, one has that, thanks to Theorem 3.1, P(F n < z) — $(z) x 
f(n) (the symbol >c indicates asymptotic equivalence) so that, for a fixed z, 
the addition of ^E(F 3 )^ 3 \z) actually increases the rate of convergence to 
zero. 



Proof of Proposition 3.3. Since assumption (a) is in order and 
E(F n ) = 0, one can deduce from Lemma 2.1, in the case 8=1, that 

l-iDF^-DL^Fn)^ 



E[F n x 



ip{n) 



1 

292(71 



■E{F 3 ). 



Assumption (b), combined with the fact that <p(n) 1 (1 — (DF n ,—DL 1 F n )^) 
has variance 1, immediately yields that there exists 5 > such that 



sup-E|F„ x ip{n 

n>l 



i-i/ 



l-{DF n ,-DL- l F 1 



1 1+5 



< OO. 



In particular, the sequence {F n x (/?(n) -1 (l — (DF n , —DL~ l F n )f ) ) :n > 1} is 
uniformly integrable. Therefore, since assumption (iii) in the statement of 
Theorem 3.1 is in order, we may deduce that, as n — > 00, 

-E(F 3 ) — -E(N!N 2 ) = -p. 



2(p(n) 



As a consequence, 
V.(n)- 1 



P(F n <z)-*(z) + ±E(F%)*®( l 



< 



P{F n <z)- Hz) 
<p(n) 



-& 3 \; 



+ 



\^ 3 \z)\ 



2<p(n 



■E(F 3 )+p 
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and the conclusion follows from Theorem 3.1. □ 

Remark 3.5. By inspection of the proof of Proposition 3.3, one sees 
that Assumption (b) in the statement may equally well be replaced by the 
following, weaker, condition: (b') the sequence 

F n x ^{n)- l {l - (DF n , —DL^Fn)^), n > 1, 

is uniformly integrable. 

3.2. Multiple integrals. The following statement specializes the content 
of the previous subsection to multiple integrals with respect to some isonor- 
mal Gaussian process X = {X(h) : h S $)}. Recall that a nonzero finite sum 
of multiple integrals of arbitrary orders is always an element of B ' 2 and, 
also, that its law admits a density with respect to Lebesgue measure [this 
implies that assumption (i) in the statement of Theorem 3.1 is automatically 
satisfied in this section]; see Shigekawa [32]. 

Proposition 3.6. Let q > 2 be a fixed integer and let the sequence F n , 
n > 1, have the form F n = I q (f n ), where, for n > 1, f n £ Sj &q . Suppose that 
E(F 2 ) = (?!||/n||i®g - ► I- Then, the quantity ip(n) appearing in formula (1.8) 
is such that 

(3.4) <p(n) 2 = E[(l-q- 1 \\DF n \\l) 2 } 

(3.5) =(l-q\\\fn\\lm) 2 

+ q 2 J2(2q - 2r)!(r - l)! 2 ( « ~_ J V ||/„® r /„|| J 9a( 
Now, suppose that, as n — > c<3 ; 

(3.6) \\f n ® r f n \\sj®2(<i-r) ~* 

for every r = 1, . . . , q — 1 and, also, 

(3.7) l-g'll/n||^ ^ a 

(p(n) 

Then, assumption (ii) in the statement of Theorem 3.1 is satisfied and a set 
of sufficient conditions, implying that assumption (iii) in the same theorem 
holds, are the following relations (3.8)-(3.9): as n^oo, 

(3.8) V{n)~ 2 \\{f n ®rfn) ®l (fn®rfn) ||jj®2(2( ff -r)-l) -> 

for every r = 1, . . . , q — 1 and every 1 = 1,..., 2{q — r) — 1 and, if q is even, 

(3.9) - M !(g/2-l)! ( ^~ } ^\(n)~ l {f n , f n ® q/2 f n )^ q — > p. 

If q is odd and (3.8) holds, then assumption (ii) in Theorem 3.1 holds with 
p = 0. 



q—r) ■ 
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Proof. Formulae (3.4)— (3.5) are a consequence of [20], Proposition 3.2. 
The fact that (3.6) implies <p(n) — ► is immediate (recall that \\f n <g> r 
/n||f,®2( 9 -r) > \\f n ®rfn\\sjm<i-r))- Again using [20], formula (3.42), one has 
that 



(3.10) 




Finally, the fact that (3.8) and (3.9) (for q even) imply that assumption (iii) 
in Theorem 3.1 is satisfied, is a consequence of representation (3.10) and 
part B of Theorem 2.6, in the case 




and c 2 = 1, by taking into account the fact that, for q even, 

E{F n G n ) = -qq\{q/2-l)\^~}^\{nr^f n J n ® q/2 f n ) mq , 

whereas E(F n G n ) = for q odd. □ 

Remark 3.7. Observe that, due to part A of Theorem 2.6, condition 
(3.6) is actually necessary and sufficient to have <p(n) — > 0. Moreover, if 
conditions (3.6)-(3.9) are satisfied, then the usual properties of finite sums 
of multiple integrals (see, e.g., [13], Chapter VI) imply that assumptions 
(a)-(b) in the statement of Proposition 3.3 are automatically met so that 
Proposition 3.6 indeed provides one-term local Edgeworth expansions. 

3.3. Second Wiener chaos. In this subsection, we focus on random vari- 
ables in the second Wiener chaos associated with an isonormal Gaussian 
process X = {X(h) : h £ fj}, that is, random variables of the type F = 12(f), 
where / G S) Q2 . To every kernel / S ft 02 , we associate two objects: (I) the 
Hilbert-Schmidt operator 

(3.11) Hf.ft^Sj; g^f®ig, 

where the contraction f ®\ g is defined according to (2.3), and (II) the 
sequence of kernels {/ ®f f :p > 1} C i} 02 , defined as follows: / (g)^ / = / 
and, for p > 2, 

(3.12) f®?f=(f® ( r 1] f)®if- 
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We write {A/j} 3 ->i to indicate the eigenvalues of Hf. Now, for p > 1, 
denote by K p (l2(f)) the pth cumulant of 12(f)- The following relation, giving 
an explicit expression for the cumulants of 12(f), is well known (see, e.g., [6] 
for a proof): one has that ki(/2(/)) = E (12(f)) = and, for p > 2, 

^(/ 2 (/))=2 p " 1 (p-l)!xTr(ifp 

(3.13) =r 1 (p-i)!x(/#" 1) /,v 

00 

= r 1 (P-i)!xEAL, 

where Tt(Hj) stands for the trace of the pth power of Hf. 

Proposition 3.8. Let F n = h(fn), n>l, be such that f n G fj 02 and 

write Kp 1 ^ = K p (F n ), p > 1. Assume that 4™^ = E(F%) — > 1 as n-t 00. 

Then, as n — > oo, -F n iV ~ ^(0, 1) i/ and on/y j/ 4 — > 0. /n i/tis 
case, we further have 



7- (n) 

(3.i4) rf K oi(^,iv) < y-^- + (4 n) - 1) 2 - 

If, in addition, we have, as n — > 00, 

«< B) - 1 



(3 15) "2 > 

( n ) /c 1 / ( n ) 1 \2 ' 

'M 76 + (4 - l) 2 

(n) (n) 

(3.16) ^ 3 — >a and — ^ ► 0, 

V«SrV6+(4 n) -i) 2 (4 n) /e+(4 n) -i) 2 ) 2 

i/ien 

P(F n < z) - $(z) a 1 , 7 2 /9 

(3.17) , ; 17 — >-^= 1-z 2 e -2 /2 osji^co. 



«r/e+(^-i) 

In particular, if a 7^ 0, i/ien i/iere exists c G (0, 1) and no > 1 snc/i that, for 
any n>iiQ, 



/ (n) 

(3.18) sup |P(F n < z) - > c J ^i- + (4 n) - l) 2 . 

2GR V 6 

Remark 3.9. 1. If E(F%) = K 2 n) = 1, then condition (3.15) becomes 
immaterial and the denominators appearing in formula (3.16) involve solely 

» 
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2. By combining (3.16) with (3.17), we have that, as n — > oo, 



(n) 

P(F n < z) - *(*) ~ ^=(1 ~ z 2 )e~ z2/2 



whenever z ^ ±1 and a ^ 0. 

PROOF of Proposition 3.8. First, since E(F n ) = 0, we have 4 = 
- 3^(^2)2. Thus, the equivalence between kJ* 5 — ► and F n ^ 
JV{Q, 1) is a direct consequence of part A of Theorem 2.6. Now, observe 
that 

\\\DF n f - 1 = 2J 2 (/„ ®i /„) + £(F 2 ) - 1 = 2/ 2 (/„ ®! /„,) + 4 n) - 1. 
In particular, 

VarQll^H 2 - l) = 8||/ n / n ||J M + (4" } - I) 2 = ^ + - I) 2 , 

where we have used (3.13) in the case p = 4 (note that (/ <g>\ /, /)^®2 = 
11/ ®i /lljsoa)- This implies that the quantity y(n) appearing in (1.7) indeed 

equals J /Q + (4™^ — l) 2 . To conclude the proof, it is sufficient to apply 
Proposition 3.6 in the case q = 2, by observing that 



1-4 W) l-2||/n||%® 



/ 4 n) /e + (4 n) -i) 2 

and, also, by using (3.13) in the casesp = 3 and p = 8, respectively, 

4^ _ 8 (fn,fn 01 /n)jj®2 



and 



4 n) /e+(4 n) -i) 2 



4 n) -2^7!x (/ ^ 7)/ ' 



(4 n) /6 + (4 n) -l) 2 ) 2 ' V(n) 4 

_ 7 . ||(/n®l/n) <8>l(/n®l/n)||%32 

4. Toeplitz quadratic functionals of continuous-time stationary processes. 

In this section, we apply our results to establish (possibly optimal) Berry- 
Esseen bounds in CLTs involving quadratic functionals of continuous-time 
stationary Gaussian processes. Our results represent a substantial refine- 
ment of the CLTs proven in the papers by Ginovyan [7] and Ginovyan and 
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Sahakyan [8], where the authors have extended to a continuous-time set- 
ting the discrete-time results of Avram [1], Fox and Taqqu [6] and Giraitis 
and Surgailis [10]. In the discrete-time case, Berry-Esseen-type bounds for 
CLTs involving special quadratic functionals of stationary Gaussian pro- 
cesses are obtained in [35], and Edgeworth expansions are studied in, for 
example, [17]. However, to our knowledge, the results which are proved in 
this section are the first (exact) Berry-Esseen bounds ever proved in the 
continuous-time case. Observe that it is not clear whether one can deduce 
bounds in continuous-time by using the discrete-time findings of [17] and 
[35]. We refer the reader to [2] and [11] (and the references therein) for 
CLTs and one-term Edgeworth expansions concerning quadratic functionals 
of general discrete-time processes. 

Let X = (X t )t£R be a centered real-valued Gaussian process with spectral 
density / : R — > R. This means that, for every u, t G R, one has 

/+oo 
e lXt f(X)dX, teR, 
-oo 

where r : R — > R is the covariance function of X. We stress that the density / 
is necessarily an even function. For T > 0, let Qt = ffi T p g(t — s)X t X s dtds, 
where 

/+oo 
e tXt g(X)dX, tern, 
-oo 

is the Fourier transform of some integrable even function g : R — > R. The 
random variable Qt is customarily called the Toeplitz quadratic functional 
of X, associated with g and T. We also set 

~ Q T - E(Q T ) - Q T 

Qt = t= and _ 



VT a(T) 

with a(T) 2 = Var(Qy). The cumulants of Qt and Qt are denoted, respec- 
tively, by 

Kj = Kj(Qt) and = Kj(QT), j > 1. 

Given T > and ip G L 1 (R), we denote by Bt(iP) the truncated Toeplitz 
operator associated with ip and T, acting on a square-integrable function u 
as follows: 

rp 

B T {ip)(u)(X)= [ u(x)$(X-x)dx, AgM, 
Jo 

where is the Fourier transform of ip. Given tp, "f G L 1 (R), we denote by 
St(V ; )-St(7) the product of the two operators Bx(ip) and Bt{^)] also, 
[Bt(i/j)Bt('j)] : ' , j > 1, is the jth power of Bt{iP)Bt{^)- The symbol Tr(U) 
indicates the trace of an operator U. 

The following statement collects some of the results proven in [7, 8]. 
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Theorem 4.1 (See [7, 8]). 1. For every j > 1, thejth cumulant of Qt is 
given by 



(T)_/0, if 3 = 1, 

T -j/^-\j - l)\Ti[B T (f)B T (g)]i , if j > 2. 



2. Assume that f G L P (M) n L 1 (IR) (p>l), that g G L 9 (R) n L 1 (R) fa>lj 
and i/iai ^ + | < | . Tnen, 

/-\-oo 
f(x) j g(x) j dx. 
-oo 

3. If l - + l<\, then 



/rn\ f+OO 

a 2 (T) = 4 T) — a 2 (oo):=16vr 3 / / 2 (x)o 2 (x) 

>oo J-OO 



and Qt ^ ~ ^(fi, 1) as T ^ oo. 



The next statement shows that one can apply Proposition 3.8 in order to 
obtain Berry-Esseen bounds for the CLT appearing in point 3 of Theorem 
4.1. Observe that, since the variance of Qt is equal to 1, by construction, to 
establish an upper bound, we need to control only the fourth cumulant of 
Qt- this will be done by using point 2 of Theorem 4.1 and by assuming that 
- + -<j. On the other hand, to prove lower bounds, one needs to have a 

precise estimate of the asymptotic behavior of the eighth cumulant of Qt- 
again in view of point 2 of Theorem 4.1, this requires that - + - < |. 

Theorem 4.2. Assume that f G W(R) n L X (R) (p > 1) and that g G 
L 9 (R) n L 1 (R) (q>l). Let<$>(z) = P(N <z), where N ~ jY '(0,1). 

1. If - + - < i, then there exists C = C(f,g) > such that, for all T > 0, 
we have 

sup\P(Q T < z) - $(z)\ < 



zgR yT 
2. If\ + \<\ and 

■+oo 



f s {x)g i {x)dx^Q, 



-oo 

then there exists c = c(f, g) > and To = 7b(/, g) such that T >T$ implies 
sup\P(Q T <z)-<S>(z)\>^=. 
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More precisely, for any z G R, we have 
Vt{P{Q t <z)-^(z)) 



(4.1) 



(1 — z )e ' . 



Proof. It is a standard result that each random variable Qt can be 
represented as a double Wiener-Ito integral with respect to X. It follows that 
the statement can be proven by means of Proposition 3.8. Now, whenever 
- + ~ < j, one easily obtains, from points 2 and 3 of Theorem 4.1, that 

(4.2) ~ t^ ^-'U ~ l)'(2vr) 2 ^ 1 J±£ f j (x)gi(x)dx 



r^oo (16vr 3 )j/ 2 (/+~ p(x)g 2 {x)dx)'Jl 2 



and the desired conclusion is then obtained by a direct application of Propo- 
sition 3.8. In particular, point 1 in the statement is immediately deduced 
from the fact that ~ + ~ < |, from relation (4.2) and the bound (3.14), with 

k^ replacing (observe that k^ = 1, by construction). On the other 
hand, point 2 is a consequence of the fact that if | + 1 < |, then (4.2) implies 
that condition (3.16) is met. The exact value of the constant appearing on 
the right-hand side of (4.1) is deduced from elementary simplifications. □ 

5. Exploding quadratic functionals of a Brownian sheet. In this section, 
we apply our results to the study of some quadratic functionals of a standard 
Brownian sheet on [0, l] d (d > 1), denoted W = {W(ii,. .-,t d ): (ti, ...,t d )£ 
[0,l] d }. We recall that W is a centered Gaussian process such that, for every 

(ti,...,t d ),(m,...,u d ) e [o,i] d , 

E[W(ti,...,t d )W(u 1 ,...,u d )]= Y[ min( Ui ,U) 

i=l,...,d 

so that, if d = l, the process W is indeed a standard Brownian motion on 
[0, 1]. It is easily proved that, for every d > 1, the Gaussian space generated 
by W can be identified with an isonormal Gaussian process of the type X = 
{X(f):f G L 2 ([0, l] d , d\)}, where dX indicates the restriction of Lebesgue 
measure on [0, l] d . It is also well known that the trajectories of W enjoy the 
following, remarkable, property: 

(5.1) / dt\ ■ ■ ■ dt d = +oo, .P-a.s. 

^[o,i]A h---t d J 

Relation (5.1) is a consequence of the scaling properties of W and of the 
well known Jeulin's lemma (see [14], Lemma 1, page 44, or [24]). In the case 
d = 1, the study of phenomena such as (5.1) arose at the end of the 1970s, 
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in connection with the theory of enlargement of nitrations (see [14, 15]); see 
also [16] for some relations with noncanonical representations of Gaussian 
processes. 

Now, denote, for every e > 0, 

and observe that B d is a centered random variable with moments of all 
orders. The CLT stated in the forthcoming proposition gives some insights 
into the "rate of explosion around zero" of the random function 




Proposition 5.1. For every d > 1 , ase^O, 
(5.2) B d := (4 log l/e)~ d/2 xB E ^JV~ ^(0, 1). 

Proposition 5.1 has been established in [27] (for the case d= 1), [5] (for 
the case d = 2) and [23] (for the case d > 2). See [27, 28] for an application of 
the CLT (5.2) (in the case d= 1) to the study of Brownian local times. See 
[5] for some applications to conditioned bivariate Gaussian processes and 
to statistical tests of independence. The next result, which is obtained by 
means of the techniques developed in this paper, gives an exact description 
(in terms of the Kolmogorov distance) of the rate of convergence of B d 
toward a Gaussian random variable. 



Proposition 5.2. For every d> 1, there exist constants < c(d) < 
C(d) < +oo and < rj(d) < 1, depending uniquely on d, such that, for every 

£>0, 

d KoX [B d £ ,N}<C{d)(\ogl/er d l 2 

and, for e < n(d), 

d Kol [B d ,N]>c(d)(logl/e)- d / 2 . 

Proof. We denote by 

Kj(d,e), j = 1,2,..., 

the sequence of cumulants of the random variable B d . We deal separately 
with the cases d = 1 and d > 2. 
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(Case d=l.) In this case, W is a standard Brownian motion on [0,1] 
so that B\ takes the form B\ = ^{fa), where 1% is the double Wiener-Ito 
integral with respect to W and 

(5.3) f £ (x,y) = (41og l/er ll2 [{x VyV e)" 1 - 1]. 

Lengthy (but standard) computations yield the following estimates: as s — > 0, 
K 2 (l,e) — ► 1, 

Kj(l,e) x ( log 7j ' J '- 3 - 

The conclusion now follows from Proposition 3.8. 

(Case d>2.) In this case, Bf has the form B^ = -^(/e ), with 

(5.4) /*(*!, . . . , x d ; yi , . . . , = (4 log l/e)"^ 2 J] [( Xj V V e)' 1 - 1]. 

By using (3.13), one sees that 

(V-\j - ly.y 1 x k^e) = [(2i-\j - l)!)" 1 x K 3 (l,e)] d , 
so the conclusion follows once again from Proposition 3.8. □ 

6. Exact asymptotics in the Breuer Major CLT. Let 5 be a fractional 
Brownian motion (fBm) with Hurst index H G (0, ^), that is, {B x : x > 0} is 
a centered Gaussian process with covariance given by 

E(B x B y ) = \(x 2H + y 2H -\x- y\ 2H ), x,y>0. 

It is well known that, for every choice of the parameter H £ (0, 5), the 
Gaussian space generated by B can be identified with an isonormal Gaussian 
process of the type X = {X(h) : h G f)}, where the real and separable Hilbert 
space Sj is defined as follows: (i) denote by $ the set of all R- valued step 
functions on M + ; (ii) define Sj as the Hilbert space obtained by closing $ 
with respect to the scalar product 

(![o,x], l[o, y ])f, = E(B x B y ). 

Such a construction implies, in particular, that B x = X(\\q^S). The reader 
is referred to, for example, [21] for more details on fBm, including crucial 
connections with fractional operators. We also define p(-) to be the covari- 
ance function associated with the stationary process x 1— ► B x+ \ — B x , that 
is, 

p(x) :=E[(B t+1 - B t )(B t+x+1 - B t+X )j 

= i(|x + l| 2 ^ + |x-l| 2 ^-2|x| 2H ), x£M. 
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Now, fix an even integer q > 2 and set 



1 



T 



Zt ■= — ~= / i?,(B u+ i - 5 U ) du, T > 0, 
a{T)\jT Jo 

where if g is the gth Hermite polynomial defined in (2.14) and where 



(r(T):= ^Vax(-^jT H q {B u+1 - B u )d^j = M j^p*{u-v)dudv. 

Observe that each Zt can be represented as a multiple Wiener-Ito integral 
of order q and also that 

/+oo 
p q (x)dx < +oo. 
-oo 

According to, for example, the main results in [3] or [9], one always has the 
following CLT: 

Law 

z T ► z~ t /K(o,i) 

T->oo 

(which also holds for odd values of q). The forthcoming Theorem 6.1 shows 
that the techniques of this paper may be used to deduce an exact asymptotic 
relation (as T — > oo) for the difference P(Zt < z) — $(z), where <&(z) = 
P(N < z) [N ~ o/K(0, 1)]. We stress that the main results of this section deal 
with the case of a generic Hermite polynomial of even order q>2, implying 
that our techniques even provide explicit results outside the framework of 
quadratic functionals, such as those analyzed in Sections 4 and 5. In what 
follows, we use the notation 

? 2 (°c) := ^f - J> - l)! 2 ( q ~ \ V (29 - 2,)! 
o- 4 (oo)^ V s - 1 / 

(6.1) x / p s (x 1 )p s (x 2 )p q - s (x 3 ) 

x p q ~ s {x2 + X3 — Xi) dxi dX2 dX3 

and 

(6.2) 7(oc) = - \ 3 ; g f / P q/2 (x)P q/2 (y)P q/ \x-y)dxdy. 



Theorem 6.1. There exists a constant C > siic/i i/tai 

(6.3) d Kol (Z T ,N) = sup\P(Z T <z)-<S>(z)\ <-^=. 

zGR V J- 
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Moreover, for any fixed z € M, we have 

(6.4) VT(P(Z T <z)-$(z)) ^(/-l)^. 

Proof. The proof is divided into three steps. 
Step 1. Let us first prove the following convergence: 

(6.5) Vf(-\\DZ T \\l - l) ^^( ,a 2 (oc)), 

\q J T-xx. 

where a 2 (oo) is given by (6.1). Note that, once (6.5) is proved to be true, 
one deduces immediately that, as T — ► oo, 

/1„ „9 \ 5 2 (oo) 
Var^-Hl^-lJ-^A 

so that (6.3) follows from Theorem 2.4. Now, to prove that (6.5) holds, start 
by using the well-known relation between Hermite polynomials and multiple 
integrals to write 

H q {B u+ i - B u ) = I q (l® u q u+1] ). 
As a consequence, we have 

q rT ^ 

DZ T = J Vi ) 1 K«+i] 

Thus, by an appropriate use of the multiplication formula (2.4), one has that 
\\DZ T \\% = y p(n - ^^-illKi])^-^?^!]) 

/ V r! ( 9 " 1 ) J 2 ,_ 2 _2r(l? 8> if 
•Ao,ri 2 4* V r y 29 2rV w [■ 



-l-r\ 



x p r+1 (u — v) dudv 



X An ^ ^(TS ® C+V^ " U ) 
I'M J 



yielding 
1 



DZ T || 2 -1 
q 
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1~ l / 1 \ 2 



I, 



'*r-*(lg£i] 8 iR+V^ " u ) dudv - 



<[0,T] 2 

We shall first prove that, for every s £ {1, . . . ,q — 1}, 



1 



/2«7-2 S (l? (9 Til 8> lf l9 7$)ff(u - v) dudv 

T J[0T] 2 [u 7 u+l] ^ [v,v+l\'r \ J 



Law o 

^(0,a s 2 (oo , 



(6.6) 
where 

a^(oo) := (2q-2s)\ / p s (xi)p s (x 2 )p q ~ s (x 3 )p q ~ s (x2 + x 3 - x x ) dx\ dx 2 dx 3 . 
Jm 3 

Fix s G {1, . . . ,q — 1}. First, observe that 

^ == Var (-L ^ /^(ig^j 8 if^V^ - «) dud,) 

= (2g ~ 2S)! / P> - v)p s (w - z)p q - s (u - w) 

J J[0,T] 4 

x /j 9-5 (i> — z) du dv dw dz 

1 — >oo 

so that (6.6) holds if and only if the following convergence takes place: 
(6.7) 

Law 

^o,i). 

We have 

Js) _ 2q-2s 



DQ 



T a s (T)VT 

x L ™ pS{u ~ ^-^-^E^ ® ^Ki ) WD 

Thus, WDQ^Wl is given by 



(2g-2s) 2 



d 2 s {T)T 7 [0 ,T]4 



p s (ui - u 2 )p s (u 3 - u 4 )p(ui - u 3 ) 
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T / 1 ®(g-s-l) x 

X i2g-2s-HJ-[ UliU1+1 ] ^ - L [ U2jU2 + l]J 

_ (2g-2s) 2 
~ 52(T)T 

X / p S («i - U 2 )p S (u 3 - U4)p(ui - U 3 ) 
J[0,T] 4 



X 



'«~ s ' 2q -2s-V 2 



t 



+ 



\t=0 

v r 1 ®(g-s-l) 
X J4g_4 s _2-2H 1 [ Ul , Ul+ l] ® ^-lua.us+l] 

..®(g-s-t) o, -i® («-*-*) \ 

^ 1 [«2,«2 + 1] 59 1 [W4,«4+1] ^ 

x /0*(n2 — ^4)^ c?u 

(2g - 2s) 2 

"WW 



x / p s (u 1 -u 2 )p s (u 3 -u 4 )p(u 1 -u 3 ) 

J[0,T] 4 



x 



\t=q— s+1 



t 



(1 2s— 1— i) 

X ^4g-4 S -2-2H 1 [ Ul , Ul+ l] 

[lt3,«3 + l] y 



x/ s {u 2 -u i )p t q+s (ui-u 3 )\dui 



Consequently, 2g -2<5 W^Qt^ llfj — 1 is given by 
2<?-2s 



x p s {ui-u 2 )p s {u 3 -u A )p{ui-u 3 ) 

J[0,T] 4 
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+ 



Xj 4g _ 4s _2HJ-[ MliU1 + i] 59 J -[w 3 ,'U3+l] 

[«2,«2+l] [«4,«4+l] ' 

x p t+l {u 2 — U/^)\ d^i • • • dli 4 

2g-2s 

p s (iii - u 2 )p s (u 3 - ua)p(ui - u 3 ) 



[0,T] 4 



\t=q-s+2 



t-1 

®(2g-2s-t) 



X i 4g -4s-2H 1 [ UliUl+ l] 



-®(2g-2s-t)<, 
[«3,U3+1] J 



x p q ~ s (u 2 - u A )p t ^i +s ~ 1 (n 1 - u 3 )J d«i • • -cb 4 . 

a fixed t such that l<i<g — s + 1, we have that 
E — ( p s (ui - u 2 )p s (u 3 - u 4 )p(ui - u 3 )p t+1 (u 2 - u 4 ) 

1 J[0,T] 4 

v r ft »(?-*- !) <> -1 ®(g-s-i) 

X i4 ? -4 S -2H 1 [ ul)U1+ l] ^ 1 [« 3 ,U3+1] 

„ ., ®(g-s+l-t) 1 ®(g-s+l-t)N j 
K)i [U2,U2+l] ^->4,«4 + l] ) aU l--- aU i 

1 



p s (ui - u 2 )p s (us - u 4 )p{ui - U 3 ) 

T- J[o,T]" 



18 

t 



X p (U2~ U4)p S (ua - U G )p S (u 7 - Ug) 

x p(u 5 - u 7 )p\u 6 - u 8 )p q ~ s ~ 1 (u 1 - U 5 ) 
x^- s - 1 (n 3 -n 7 )p g - s+1 -*(n 2 -u 6 ) 
x p 9 ~ s+1_ *(u 4 — us) du\ ■ ■ ■ du$ 

— [ p s {xi)p s (x 2 )p(x 3 )p\x 2 + x 3 - ^i)p s (x 4 )p s (x 5 )p(x 6 ) 



X p t (x 5 + X6-X 4 )p q S 1 (X 7 ) 

x p q - s - 1 (x G + x 7 - x 3 )p q ~ s+1 ~ t (x 4 + x 7 - Xl) 
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X /9 9 ~ S+1 ~ < (X5 + Xq + Xj — X2 — Xs) &X\ • ■ • dxj 

tends to zero as T — > oo. Similarly, we can prove, for a fixed t such that 
q-s + 2<t<2q-2s-l, that 

E — / p s ( Ul - u 2 )p s (u 3 - u 4 )p*" 9+s (ni - uz)p q ~ 8 {u 2 - u 4 ) 

X /4,-4 fl -2t(l [uiiUl+ i] ® 1 [«; j „ 3+1 ] )dUl---d«4 

tends to zero as T — > oo. Thanks to the main result in [22], the last relation 
implies that, for each s, the convergence (6.7) holds and, therefore, (6.6) 
is verified. Finally, by combining (6.6) with the results in [25] and [26], we 
obtain (6.5). Indeed, by using the orthogonality and isometric properties of 
multiple stochastic integrals, we can write 



Var(VrQ||J3Zr||J-l)) 



" E(s-1)< 2 (!-!) (2*- 2-)! 



'[0,T] 2 

[0,n a 1 ' + 1 1 ' + 1 /i3®(2 9 -2a) 



„2 9 -1 / 1 \ 4 

9 »-i) !2 (":l) < 2 «- 2s ) ! 

p s (u — v)p s (w — z) 



'[0,T] 4 

x p q ~ s (u — w)p q ~ s (v — z) dudvdwdz, 
from which we easily deduce that Var( v / T( i \\DZ T \\\ - 1)) — > ? 2 (oo). 

Step 2. Let us prove the following convergence: 
(6.8) ^ t ,VtQ||DZ t ||2 _i^Jf^(t/,y) 

with (U,V) a centered Gaussian vector such that E(U 2 ) = 1, 

E(V 2 ) = a 2 {oo) 
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and 

E{UV) = -7(00) 

= g!( o / 3 ) / (g f )2 / P q/ \x)P q/2 {y)P q/ \x-y)dxdy. 

la 6 (00) J R 2 

Observe that we already know that I ^ [/ and also that (6.5) is verified. 
Note, also, that we have proven (6.5) by first decomposing VT(q~ 1 \\DZT\\j d — 
1) into a finite sum of multiple integrals and then by showing that each mul- 
tiple integral satisfies an appropriate CLT. As a consequence, according to 
part B of Theorem 2.6 [with G n replaced by \/T(q~ 1 ||DZt||« — 1)], it is 
sufficient to show the following convergence: 

e(z t xVt(-\\dz t \\1-i 

(6 ' 9) 1 2 , , 2 



T^oo 2<7 3 (oo) 

By the orthogonality of multiple stochastic integrals, we can write 

e(z t x Vt(-\\dz t \\%-i 



/ 1 \ 2 



z x / n — 1 \ 2 ^ 

(Pl 2 {u- v)p q l 2 {u-w) 



T] 3 

x p q l 2 (w — v) du dv dw 



/ — 1 \ 2 

(I " i) ! y L 1 J f R2 P q/2 {x)p q/2 {y)p ql2 (x - y) dx dy 



T^oo a 3 (oo) \2 

= g !( g /2)!(J 2 ) r p^ {x -y) dxd y. 

2cr d (oo) Jm.2 

Step 3. Step 1 and convergence (6.8) imply that, as T — > 00, 
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where <p(T) = Var(l - q-^DZ^W - 1) and 
/ l-l/g\\DZ T \\l \ Law 

As a consequence, one can apply Theorem 3.1 in the case p= (the 
remaining assumptions are easily verified), yielding that 

V (T)-\P(Z T <z)-^(z))-^ J^^-l)L^. 

T->oo 6a(oo) a/27t 

from which the conclusion follows. □ 
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